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ABSTRACT 



We present a re-analysis of the Geneva-Copenhagen survey, which benefits from the infrared flux method to improve the accuracy 
of the derived stellar effective temperatures and uses the latter to build a consistent and improved metallicity scale. Metallicities are 
calibrated on high-resolution spectroscopy and checked against four open clusters and a moving group, showing excellent consistency. 
The new temperature and metallicity scales provide a better match to theoretical isochrones, which are used for a Bayesian analysis 
of stellar ages. With respect to previous analyses, our stars are on average 100 K hotter and 0.1 dex more metal rich, which shift the 
peak of the metallicity distribution function around the solar value. From Stromgren photometry we are able to derive for the first 
time a proxy for [ff/Fe] abundances, which enables us to perform a tentative dissection of the chemical thin and thick disc. We find 
evidence for the latter being composed of an old, mildly but systematically alpha-enhanced population that extends to super solar 
metallicities, in agreement with spectroscopic studies. Our revision offers the largest existing kinematically unbiased sample of the 
solar neighbourhood that contains full information on kinematics, metallicities, and ages and thus provides better constraints on the 
physical processes relevant in the build-up of the Milky Way disc, enabling a better understanding of the Sun in a Galactic context. 

Key words. Stars: abundances - Stars: fundamental parameters - (Stars:) Hertzsprung-Russell and C-M diagrams - Stars: kinematics 
and dynamics - Galaxy: disk - (Galaxy:) solar neighborhood 



1. Introduction 

Late-type dwarf stars are long-lived objects and can be regarded 
as snapshots of the stellar populations that are formed at differ- 
ent times and places over the history of our Galaxy. Not only 
their kinematics carry residual information on their dynamical 
histories, but their atmospheres retain a fossil record of the com- 
position of elements in the interstellar medium at the time and 
place of their formation. Therefore, F, G, and -to a lesser extent- 
K dwarfs have been traditionally used to study various aspects of 
the chemical evolution of the Milky Way. 

The region in the Milky Way for which this task can be 
most easily achieved is the solar neighbourhood; starting from 
pioneering works using spectra or ultraviolet and colour ex- 
cess to est imate the metal abundance of stars in a Galactic con- 
text ( e.g.. IWallersteinI 119621; Ivan den Ber"ghlll962l; lEggen et al.l 
119621; ISchmidtlll963h . this endeavour has continued over the 
years with steadily improving spectroscopic and photometric 
studie s. The latter (e.g., Twarog 1980; Olsen 1983; Stromgren 
ll987L iNordstrom eT al.1 l2004t iHavwoodl 12008 ) comprise large 
catalogues, but have to pay for this by being only able to derive 
one single parameter for metallicity, and no detailed elemental 
abundances. On the other hand, spectroscopic studies are still 
limited to small samples of a few hund red or about a thousand 
stars at most. While some sttidies ( e.g., lEdvardsson et al.ll 19931: 
iFavata et al.lll997t lFuhrmannll2008l) rely on kinematically unbi- 



ased samples, many inve stigations (e.g., Feltzing & G ustafssonl 
'l9 98tlBensbv et al. 20 03l;lReddv et alj2006l;lBensbv et al. 2003 



iRamfrez et al. l2007i; iSoubiran et al.l 2008') make use of sophis- 
ticated kinematic selections to achieve significant numbers of 
members belonging to different subpopulations in their sample. 
Even though the abundance trends in these studies are better 
traced thanks to this strategy, a quantitative interpretation can 
be more difficult. 

Galactic chemo-dynamic stu dies are now entering a new 
realm with curre nt (e.g., RAVE ISteinmetz et al.l l2006l SDSS 
Ivezic et al. 2008) and forthcoming (e.g., SkyMapper, APOGEE, 
HERMES, LSST, Gaia) large photometric, spectroscopic and as- 
trometric surveys targeting different and fainter components of 
the Galaxy. These tremendous observational efforts, however, 
must be supported by equal investments to minimize the eiTors 
that plague the determination of stellar parameters. The most im- 
portant parameter is the effective temperature (reft): its deter- 
mination has implications for the derived abundances, for sur- 
face gravities and for the inferred ages, masses, and distances 
via isochrone fitting. If we aim to deconstruct the formation and 
evolution of the Milky Way in a star-by-star fashion, it is funda- 
mental to have full control over all potential sources of errors. 

The prefeiTed stellar Teff scale has been a long debated issue, 
with various scales differing systematically by 100 K or more. 
Though this is still true in many areas of the HR diagram, re- 
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cent d ata on solar twins (iMelendez et alJ l2009t iRanurez et alJ 
12009^. new data and analyses of interferometric angular diam- 
eters (Bovaiian et al. 2010; Chiavassa et al. 2010) and improved 
HST absolute spectrophotometry (Bohlin 2007) have allowed to 
pin down the source of these discrepancies via the infrared flux 
method (IRFM). This gives a good base for the zeropoint of the 
temperature scale of dwarfs and subgiants, which has now an 
uncertainty of the order of only 20 K (Casagrande et al. 201^. 
For solar-type stars, the new IRFM scale supports eifective tem- 
peratu res approximately 100 K hotter than those of lAlonso et alJ 
(Il996l) . which has been the de facto choice in many studies un- 
til now. Such a shift on the zeropoint has an immediate conse- 
quence on the a bundances and ages der ived for nearby, solar- like 
stars (see also IMelendez et al.l 2010bh and therefore for inter- 
preting the most basic constraints on Galactic chemical evolu- 
tion models, namely the metallicity distribution function and the 
age-metallicity relation. The HR diagram constructed using our 
newly derived T^s scale matches very w ell that predicted by stel- 
lar models for evolved F and G dwarfs (IVandenBerg etal.ll2010i 
iBrasseur et al. l l20T0l) . thus avoiding the introduction of any ad 
hoc shift s to the Tefr scale as was the case in some previous stud- 
ies re.g.. lNordstr6m et al.ll200 4). 

The purpose of the present work is to carry out a revision of 
the astrophysical pararn eters in the Gene va-Copenhag en Survey 
(lNordstr5m etal.ll2004 iHolmberg et aH |2007, 2009) with the 
new e ffective temperature scale presented in Casagrande et ak] 
( 1201 Ol) as a starting point to derive new metallicities and ages. 
We improve not only on the accuracy, i.e. reduce zeropoint sys- 
tematics, but also the precision by reducing internal errors stem- 
ming from photometric transformations, resulting in highly ho- 
mogeneous astrophysical parameters. These improvements turn 
out to be crucial to provide more stringent observational con- 
straints on Galactic chemical evolution theories and hence on the 
history of the Milky Way. In fact, a knowledge of the metallic- 
ity distribution together with Galactic abundance gradients can 
improve our understanding of the impact a nd shape of the stel- 
lar migr ation process in the Galactic disc (iSchonrich & BinnevI 
I2009atbl) . Because models including radial migration relax the 
classical tight correlation between age and metallicity, this re- 
lation becomes effectively an additional constraint independent 
from the metallicity distribution. 

As we will demonstrate, an estimate of [a/Fe] for most of 
the stars is also obtained here for the first time from Stromgren 
indices. Having an indication of [a/Fe] allows for a tentative dis- 
section of the chemical thin and thick disc. These estimates are 
far less accurate than those obtained by high-resolution spec- 
troscopy, yet this sample exceeds the largest spectroscopic stud- 
ies available so far by more than an order of magnitude and it is 
not biased by any kinematic selection. 

The paper is organized as follows. We present the sample and 
the determination of new effective temperatures and metallici- 
ties in Section 121 Correspondingly, new ages and masses for the 
stars are derived in Section |3] In Section |4] we use this informa- 
tion for studying the metallicity distribution function in the solar 
neighbourhood and briefly discuss a possible signature of the 
Galactic bar. The age-dispersion relation is discussed in Section 
|5] while Sectionals devoted to a better understanding of the disc 
and its metallicity gradient. We finally present our conclusions 
in Section 121 



2. Determination of astrophysical parameters 

The Geneva-Copenhagen Survey (GCS) is the most compre- 
hensive catalogue of late-type solar neighbourhood stars, pro- 
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Fig. 1. Distributions in the velocity planes of the two subsamples 
defined in this work. 



viding kinematics and Galactic orbits for a magnitude-limited 
and kinematically unbiased sample of 16682 of FG(K) dwarfs 
brighter than V ~ 8.3. Some 63000 radial velocity measure- 
ments were used to assemble the catalogue, which, comple- 
mented with Tycho2 proper motions and Hipparcos parallaxes, 
also provides binaiity indication and distances. Because the se- 
lection of stars into the final catalogue was purely based on 
colour and magnitude cuts, the survey provides (apart from ef- 
fects by the photometric selection) a kinematically unbiased cen- 
sus of the solar neighbourhood. While we do not have access to 
the original sampl e selection performed in assembling the cata- 
logue, we refer to Nor dstrom et al 1 (12004 for a comprehensive 
discussion on the adopted selection criteria and relative com- 
pleteness. Homogeneous Stromgren photometry was used to de- 
rive Teff and [F e/H] for nearly a ll stars in the survey. The orig- 
inal catalogue dNordstrom et al.l [2004, GCSI) has undergone a 
number of re visions to improve th e temperature and metallicity 
calibrations ( Holmberg et al.l 12007. GCSII) an d to account for 
the new reduction of the Hipparcos parallaxes (iHolmberg et al.l 
2009, GCSIII). 

However, recent work has show n that the temperature 
scale ad opted in G CSI-III is too cold dCasagrande et al J 120101 : 
Melendez et al.ll20l0bi) . This has far-reaching impUcations: hot- 
ter temperatures imply higher spectroscopic metallicities and - 
when relying on stellar isochrones- lower age estimates. 

Because we use photometry to derive astrophysical parame- 
ters, it is crucial to clean the sample from binaries, variable stars 
and/or less certain measurements. A description of our selection 
leading to stars with the best photometry {irfm sample) with re- 
spect to the remaining ones {clbr sample) is given in Section lTTI 
where we also briefly present our implementation of the IRFM 
and the new effective temperatures derived for the entire GCS 
catalogue. Notice that the distinction between the two samples is 
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Fig. 2. Panel a) and b): AT^ff (ours minus GCSII) as a function of our Teff (upper) and [Fe/H] (lower) for stars in the irfm sample. 
Panels c) and d): same as before, but for all remaining GCS stars in the clbr sample (see discussion in Section lTTT l. Contour levels 
are computed on abscissa intervals of 25 K and 0.025 dex, respectively, to equally represent regions with fewer stars. Dashed and 
dotted lines indicate the mean difference and standard deviation of the entire sample; because they are symmetric and are dominated 
by regions with the highest overdensity of stars, the dashed and dotted lines are in some cases offset from the local Icr contour 
levels. 



based exclusively on the photometric criteria applied and there- 
fore does not introduce any apparent kinematic bias (Fig. [T). The 
corresponding new metallicity scale and ages are then discussed 
in Section lZ2l and[3l respectively. 



2.1. A new effective temperature scale 

The effective temperature s in the GCSI were derived using the 
Stromgren calibration of l Alonso et al.l (Il996l) . which however 
lacked a sufficient number of stars bluer (i.e. hotter) than (b-y) ~ 
0.3 (T^eff ~ 6500 K). In GCSII this problem was tackled by de- 
riving a new (b-y) vs. Teff calibration, where effective tempera- 
tu res for all stars were first obtained using the (V-K) calibration 
of Idi Benedettol (11998.) and then T^g were re-derived by apply- 
ing this new (b - y) cahbration to the full catalogue. However, 
the calibration of I di Benedettol (Il998h is defined in Johnson K, 
enforcing a colour transformation from the 2MASS Ks used in 
GCSII. Because the standards of the Johnson system are all satu- 
rated in 2MASS, this renders the transformation between the two 
systems less precise (ICarpenterl l200 1 ) . In addition, the metallic- 
ity effect is largely reduced but not zero even in (V - K), and the 
calibration of di Benedetto (1998) does not account for this ef- 
fect. The zeropoint of the ldi Bene detto ( 1998) t emperature scale 
is roughly intermediate between that of the Alonso et al.l(ll996l) 



scale and our own, which is some 50 K hotter (see below and the 
comparison in lCasagrande et ani2006h . 

2.1.1. Reddening 

When deriving Teff from photometry it is crucial to correct for 
reddening, if present. Fortunately, given the solar neighbourhood 
nature of the sample used here, most of the stars are unaffected 
by this problem. Reddening estimates derived from Stromgren 
photometry are known to be generally reliable (in this case with a 
stated p recision crE(b-v) - 0.009 m ag, .Holmberg et al. 2007, but 
see also lKarata^ & Schusteill20Tol for a recent revision), and we 
adopted a procedure similar to GCS for all starsQi.e. a reddening 
correction is applied for stars with E{b - y) greater than 0.01 



' In the irfm sample the col our excess has been sc aled according to 
intrinsic colour of the star ("see lCasagrande et al.ll2010i) from which the 
following mean reddening relations were computed and used for the 
clbr sample: E{Bt - Vj) = 1.32 E(b - y), E(Vt - J) = 3.18 E(b - y), 
E(Vt-H) = 3.66 E(b-y),E(VT-Ks) = 3.93 E(b-y) and Rbj =4.23, 
Rvj = 3.24, Rj = 0.86, Rh = 0.50, Rk, = 0.30, where R^ = -0^^, 
A(^) is the extinction in a given ( band. For the Stromgren indices 
E(mi) = -0.30E(b - y) an d E(ci) = 0.20 E(b - y) were adopted from 
ICrawford & BamesI (Il970l) . 
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(panel b). Contour levels and lines as in Fig. |2] 



mag and farther away than 40 pc, otherwise a value of zero is 
assumed. 

Thus, only about one quarter of the stars in the GCS cata- 
logue need to be corrected for reddening, the median E{b - y) 
being 0.02 mag, as one would expect given t he nearby nature 
of the sample (see also iHolmberg et al.ll2007l for a plot of the 
colour excess in different distance intervals). Note that the effect 
of colour excess on Teff derived via IRFM amounts to ~ 50 K for 
every 0.01 mag (see ICasagrande et al. 2010, for further details) 
whereas in the case of colour-temperature caUbrations its effect 
can be directly estimated. 

Although the temperature and metallicity scales are un- 
changed between GCSII and III, we noticed important differ- 
ences between the two catalogues. These differences show a cor- 
relation with the adopted E{b - y), reaching several hundred 
K in Teff and almost 1 dex in [Fe/H] for stars with the high- 
est colour excesses. This suggests that stars in GCSIII have not 
been corrected for reddening and because of this we will only 
use the kinematic data from GCSIII and the stellar parameters 
from GCSII when making comparisons. 

2.1.2. irfrn sample 

The IRFM implementation described in'Casagr ande et aO (1201 Oh 
not only improves the accuracy of the zeropoint of the derived 
stellar parameters, but also their precision by employing Tycho2 
BtVt and 2MASS JHKs photometry to simultaneously recover 



the bolometric flux -;Fso/(Earth)- and the effective tempera- 
ture of each star It is well suited to be applied to the Geneva- 
Copenhagen catalogue directly, avoiding the use of colour cali- 
brations as well as transformations among different photometric 
systems. 

Given its nature, it is crucial to have good photometry in 
all bands for the stars we apply the IRFM to. From the GCSII 
we exclude stars flagged as variable or having multiple com- 
ponents. \ye r etrieve Tycho2 BtVt magnitudes for all targets 
(lH0g et ani200 0) and additionally cross-check and discard those 
classified as variable or non-single in Hipparcos. The faintest 
stars might have uncertain photometry in Tycho2, whereas the 
brightest can be saturated in 2MASS: when applying the IRFM 
we consider only stars with photometric errors ctbj +crvj < 0.10 
and "j_"-i-"h_"-i-"k_msigcom"< 0.10 all having quality flag "A" in 
2MASSB Stars having T^ff < 5000 K emit considerable amount 
of flux in the red. The computation of the bolometric flux (and 
thus reff) is inaccurate if one use s only Tycho2 and 2MASS pho- 
tometry (ICasagrande et al.l|2010h and therefore we exclude stars 
cooler than this limit. This cut concerns only a minor part (326 
stars out of 16682) of the sample in the GCS. 

In our implementation of the IRFM an iterative procedure 
is adopted to cope with the mildly model-dependent nature of 
the bolometric correction: given an initial estimate of Teff, we 
interpolate over a grid of synthetic stellar fluxes at the appropri- 
ate [Fe/H] (as determined in Section |2!2] | and \ogg of each star, 
until convergence in Tetf is reached within 1 K. 

For all stars, log g is determined from the fundamental rela- 
tion 

log^=log^+41og^-logil, (1) 

where L is the bolometric luminosity^ and M is the mass of the 
star, obtained by interpolating over isochrones. Notice that in 
Eq. ([T]i mass plays only a secondary role: varying it by 10% 
changes \ogg by 0.04 dex. We used the masses reported in 
GCSII as a starting value and the BASTI mass expectation val- 
ues (cf. appendix) for a second iteration. Even variations as 
large as +0.5 dex in surface grav ity change the T^g obtained vi a 
IRFM by only a few tens of a K (Casagrande et alj|2006ir201(]h . 
thus having negligible impact. The bolometric luminosity L is 
computed from 'FB„/(Earth) using the new Hipparcos parallaxes 
(van Leeuwen 2007), and an iterative procedure was adopted to 
converge in log g using at each step the corresponding effective 
temperature and luminosity obtained from the IRFM. Although 
in the GCSI a photometric selection was made to cut out giant 
stars, there is a handful of them left. We exclude those labelled 
as suspected giants in the GCS and restrict the irfin sample to 
iogg > 3.0. Altogether, we are left with a sample of 6670 stars 
that satisfy all of the above criteria on photometric quality, non- 
binarity, and surface gravity. To these we can apply the IRFM. A 
MonteCarlo simulation using the measured observational errors 
{(TB-r, (Tvj, "j-", "h_", "k_msigcom" and cr[pe/H]) was used to es- 
timate the random error in the resulting T^s and ;FB„/(Earth) of 
each star, to which the systematic uncertainty arising fr om the 
adopted ab solute calibration was added (see .Casagrande et al.l 
l2006ll20T0h . 



^ i.e. with best photometric detection 'http:// www.ipac.cal tech.edii/l 
2mass/releases/allsky/doc/secl_6b.html#phquaI 
^ In this work we take Lq = 3.842 x lO'^'ergs-' jBahcall et all2006h . 
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2.1 .3. clbr sample 

For all remaining stars in the GCSII effective temperatures and 
bolometric fluxes were computed using the colour calibrations 
in {b - y), {Bt - Vj), ( Vt - J), (Vt - H) and (Vt - Ks) from 
ICasagrande et all (1201 Ol) . which extend also below 5000 K. We 
only took into consideration photometry with cr^^ < 0.05, o-v-^ < 
0.05, "j_msigcom"< 0.04, "h_msigcom"< 0.04, "k_msigcom"< 
0.04 (which in the colour-temperature relations imply formal 
uncertainties similar to those of the stars analysed using the 
IRFM). We computed the average Teff and ;FBo/(Earth) if more 
than one index was used and applied a 3cr clipping if more than 
two indices were present. In the latter case, the standard devia- 
tion was used to estimate the error in the derived T^ff. Notice that 
these calibrations (as any available in literature) do not include 
an explicit log g dependence. However, because surface gravities 
of dwarfs and subgiants decrease when moving to hotter Teff , an 
intrinsic dependence on su ch a term is likely to be built into them 
(see also the discussion in lCasap"ande et al . 2010). Our colour- 
temperature calibrations indeed perform extremely well along 
most of th e CMP morphology defi ned by F and G dwarfs and 
subgiants (IVandenBerg et aU20l (f). 

Figure |2] shows the comparison between T^ff derived in 
Section 12.1.21 and 12.1.31 and those in GCSII. A mean differ- 
ence of about lOOK appears, and there are trends at the high- 
est and lowest T^g, as well as at the lowest metallicities. The 
latter trend could a rise from the absenc e of an explicit metal- 
licity dependence in Idj Benedettol iL99.8|) or from the limitation 
of the standard functional form used in literature when fitting 
effe ctive temperatures and metallicities as function of (b - y) 
(see ICas agrande et al. 2010). We note that the IRFM depends 
only marginally on the assumed metallicity, and we verified for 
the GCS stars that changing [Fe/H] by ±0.2 dex affects T^ff by 
about 20 K at most. The impact can indeed be larger when one 
uses colour-reff relations that involve optical bands. 

2.2. A new Stromgren metallicity scale 

The uvby photometricQ system dStromgrenll 1 963h is well suited 
for the determination of basic stellar atmospheric parameters 
through the colour indices (b - y), m\ - (v - b) - (b - y) and 
c\ - (m - v) - (v - b). The m\ index is designed to measure 
the depression owing to metal lines around 4100 A (v band), 
and h ence is suitab le to infer the metallicity in a variety of stars 
(e.g.. lBessellll2005l and references therein). The c\ index is de- 
signed to evaluate the Balmer discontinuity, which is a temper- 
ature indicator for B- and A-type stars and a surface gravity 
indicator for late-type stars, though for stars comparable to or 
cooler than the Sun it also carries meta llicity information (e.g ., 
iTwarog et al.ll2002t lOnehag et al.ll2009t iMelendez et ani2010bi) . 
Several calibrations exist in the literature that link Stromgren 
colours to astrophysical parameters, following either theoretical 
(i.e. based on model atmospheres) or empirical approaches (see 
IOnehagetaT]l2009l and lArnadottir et all I20T0I respectivelv. for 
recent reviews). 

Throughout the paper, we talk of metallicity both in terms 
of iron abundance [Fe/H] and overall metal content [M/H], if 
a clear distinction is not needed. Notice though that Eq. ^ 
and (O given later in this section are calibrated against spec- 



In the following, we will refer to {b — y), mi and ci with the implicit 
understanding that they were dereddened if there was any colour excess. 
In the same manner, absolute magnitudes were corrected as well when 
necessary. 



troscopic measurements of [Fe/H] and thus are strictly referring 
to iron abundance. The overall metal content -always indicated 
by [M/H] in this work- was obtained using the same functional 
forms, but accounting for [a/Fe] in the calibration sample. Later 
in this section we develop a new estimator for the a-element 
content in most of the GCS stars. 

The metallicity calibration adopted in GCSII patches the red 
(b-y > 0.46) and blue (b-y < 0.30) calibrations derived in the 
GCSI with a new calibration containing all possible combina- 
tions of (b—y), m\ and ci to third orderfor0.30 < {b-y) < 0.46. 
Those three calibrations are built by linking Stromgren indices 
to spectroscopic metallicities gathered from a large number of 
studies, affecting the homogeneity of the results. In addition, the 
calibration in the 0.30 < {b -y) < 0.46 range is based on spec- 
troscopic studies with a T^s scale broadly consistent with that 
adopted in the GCSII, i.e. cooler than the one used in this study, 
implying an offset in the zeropoint of the metallicity scale. The 
adopted Teff scale is in fact the main driver in setting the zero- 
point of the metallicity scale. 

Over the past few years an increasing number of high- 
resolution and high signal-to-noise spectroscopic investigations 
have targeted hundreds of stars in the solar neighbourhood. This 
allows us to build a large and homogeneous spectroscopic cat- 
alogue, which we use to derive a new metallicity calibration. 
To this purpose we have taken only three large surveys, namely 
Valenti & Fischer (2005, V05), Sousa et a l. (2008. S08) and 
Bensby et al. (201 1 in pre p., Bl 1, which inclu des over 600 stars 
in addition to 102 from iBensbv et al.l l2003l and iBensbv et al.l 
2005). Apart from spectroscopicafly determined Teff and [Fe/H], 
all these surv eys provide a abundances: Si and Ti in the case 
of Valen ti & Fis cher (2005), and Mg, Si, Ca, Ti for the other 
two studies (for the Sousa et al. 200 8 sample the abun dances 
are given in the companion paper of iNeves et al.ll2009l N09). 
They are all very consistent, with mean differences (all in the 
sense B11-V05 and B11-S08 for 142 and 85 stars in common, 
respectively) of A[Fe/H] = 0.034 + 0.004 (cr = 0.050 dex) and 
0.047 ± 0.005 (o- = 0.046 dex) and A[a'/Fe] = 0.06 + 0.01 
(o- = 0.16 dex) and 0.01+0.02(0- = 0.20 dex). These differences 
are small and consistent with the scatter; we also made an at- 
tempt to homogenise all stars on a common scale (Bl 1) by fitting 
the differences with respect to B 1 1 as linear or parabolic function 
of [Fe/H], \ogg and T^s but this approach only had a minor ef- 
fect on the overall metallicity calibration. A comparison with the 
homogenised spectroscopic catalogue of Arnadottir et al. (201^ 
confirms this conclusion (see below). 

Our final sample contains 1522 stars, all having Stromgren 
colours, [Fe/H] and [a/Fe]. If a star was found in more than one 
study, we chose the [Fe/H] and [c/Fe] from the one that had 
Teff closest to our estimate. The mean difference between pho- 
tometric and spectroscopic T^g is 13 + 95 K. We also applied 
a 3cr clipping to remove the major outliers, and obtained a fi- 
nal calibration sample of 1498 stars (Aretf = 14+83 K), half 
of which are in the irfm sample. Fig. |3] compares our effective 
temperatures with those of the three spectroscopic studies. The 
systemati c offset between older photo metric and spectroscopic 
Teff (e.g., iRanurez & Me lendez' 2004) is now clearly removed 
thanks to our new IRFM implementation. There are no signif- 
icant trends as function of effective temperature, except for the 
very few stars below ~ 5000 K where spectroscopic estimates 
have the tendency to r eturn hotter Jeff th an photometric ones (see 
also the discussion in lSousa et al.ll2008l) . When plotting ISJ^s as 
function of [Fe/H], the metal-poor stars are on average well re- 
produced despite an increasing scatter There is a minor trend in 
the range -0.5 < [Fe/H] < 0.5 dex, with Icr contour going from 
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Fig. 4. Left panel: [Fe/H] versus mi for our 1498 calibrating stars in different (b-y) ranges represented by different symbols (median 
values as indicated in the labels). The size of the symbols increases with higher values of ci. Dotted and continuous lines represent 
Eq. ^ at these median values for the 10th and 90th percentiles of the ci distribution. Central panel: same as left panel, but for 
Eq. (O, which applies only to cool stars. Right panel: theoretical [Fe/H] versus nii relation when using synthetic colours from the 
"MARCS-standard" library and log^ instead of ci (see text for discussion). 



+50 K to -50 K: this could potentially introduce a mild system- 
atic bias (as well as affect the width of the metallicity distribu- 
tion function) of the order of +0.05 dex throughout this range, 
though for a single star this is below the accuracy of our calibra- 
tion (see below) and spectroscopic measurements are themselves 
not immune from deficiencies. On average there is no significant 
zeropoint offset or trend. 

In the literature various approaches have been used to 
calibrate Stromgren photometry to derive metallicities, ei- 
ther based on how much the colour indices mi and ci 
differ from a gi ven standard relation (usually derived for 
the H yades, e.g., lOlsenI n"984l : iHavwoodl l2002t Ivan LeeuwenI 
l2009h or using direct combinations of the Stromgren in- 
dices mi,ci and (b - y). This is the choice made in 
most of the r e cent works (e.g ., Schuster & Nissen 19 891 
Haywood 2002; Nordstrom et al. 2004; Ram irez & MelendezI 
2005: ^olmberg et al. 2007: Twarog et al. 2007h . We adopt the 



latter approach, but we are aware that even though our caU- 
brating sample includes a large number of stars, some regions 
of the [Fe/H], Tgif and \ogg space are less well sampled than 
others (see also Fig. [Hi. To limit possible biases, we checked 
our findings against synthetic colours. Despite the inaccura- 
cies that might still p lague synthetic Stromgren colours (e.g., 
iMelendez et al."2010b'). in many cases they can be used at least 
to provide guidance on general trends (Qnehag et al. 2009). For 
this work, synthetic indices were co mputed for the full gri d of 
"MARCS-standard'Q model spectra (iGustafsson et al.ll2008l) us- 
ing the zeropoints and filter transmission curves described in 
IMelendez et al. ( 2010b). Note that the purpose of using synthetic 
colours is for verification only, and they do not enter into our cal- 
ibrations, which remain fully empirical. 



^ http://marcs.astro.uu.se where standard refers to the chemical com- 
position, i.e. [cy/Fe] = 0.0 for [Fe/H] > 0.0, a linear increase of [a/Fe] 
from O.I at [Fe/H] = -0.25 to [a/Fe] = 0.4 at [Fe/H] = -1.0 and 
[a/Fe] = 0.4 for [Fe/H] < -1.0. 



Fig. m shows the sensitivity of mi to [Fe/H] for our 1498 
calibrating stars, in different {b - y) (basically Tefj) ranges. The 
asymptotic behaviour towards the most metal-poor stars reflects 
the decreasing sensitivity of mi in this regi me and it can be well 
represented by a logarithmic term (Sch uster & Nis sen 1989}. 
Therefore, rather than including all possible combination of in- 
dices in some high-order polynomial, we started with a simple 
form of the kind log(mi) + a mj and introduced mixed terms to 
allow for a change of slope with (b - y) and ci, where the ra- 
tio between the logarithmic and cubic terms, a, was optimized 
by treating it as a free parameter in the fitting process. This ac- 
counts for the first six terms in the following equation 

[Fe/H] = 3.927 log(OTi) - 14.459 mi^ - 5.394(^-3;) log(mi) 
+ 36.069 (b-y)mi^ + 3.537 cj log(mi) 
- 3.500 ot]ci + 11.034 (b-y)- 22.180 (b-yf 



+ 10.684 ci -6.759 cf- 1.548 



(2) 



where the additional terms that have a linear and quadratic de- 
pendence on (b - y) and ci were introduced after verifying that 
they improved the residuals. We also checked that the inclusion 
of terms of higher order did not lead to any further gain. Equation 
(|2]l applies to stars in the following ranges: 0.23 < (b-y) < 0.63, 
0.05 < mi < 0.68 and 0.13 < ci < 0.60 with a standard de- 
viation of 0.10 dex. We remark that for stars with [Fe/H] < 
-2 Stromgren indices effectively lose sensitivity to metallicity. 
We verified this u sing an additiona l sample of 26 metal-poor 
dwarfs taken from lCasagrande et al.l (|2010) and M elendez et al.l 
(I2010al) . Those stars, afl in the range -3.3 < [Fe/H] < -2.0, did 
not show any significant dependence on metallicity and were 
therefore not used in the fitting process, which was limited to the 
1498 stars shown in Fig.|4] 



* This appears not to be the case for very metal-poor giants above 
the horizontal b ranch where in fact mi follows [Fe/H] tightly (e.g. 
I Aden et aI.ll201lL Aden et al. to be submitted). 
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Fig. 5. Panel a): spectroscopic versus photometric metallicities 
obtained using the caUbration presented in Section IZ2l for 1498 
stars. Filled circles are for stars having {b - y) < 0.43, open cir- 
cles for redder (i.e. cooler) stars. Panel b): same as above, but 
showing residuals (ours minus spectroscopic). Dotted lines are 
Icr scatter and boxes the median values in non overlapping in- 
tervals of 0.2 dex. Panel c): same as in the previous one, but as 
a function of Teff, with boxes computed in non overlapping in- 
tervals of 100 K. The overall zeropoint offset is yu = -0.003 dex 
and cr - 0.097 dex. Panel d): distribution of the residuals of our 
calibration against spectroscopy with a Gaussian of width cr and 
centred at overplotted. Panel e): standard deviation associated 
to each square computed in panel b) with error bars being the 
standard deviation of the mean. 



While the hottest stars in Fig. |4] display a remarkably 
tight correlation with [Fe/ H], for decreasing T ^ s also ci corre- 
lates w ell with metallicity (Twarog et al. 2002; 'Melendez et al.' 
1201 Obi) . In fact, different metallicity calibratio ns are often given 
for F and GK dwarfs separately (e.g iSchuster & Nissenlll989l: 
iNordstrom et al.ll2004l) . While Eq. ([2) applies also to cool stars 



((b -y) > 0.43, i.e. Tes < 5600 K), for those we found an addi- 
tional function of the kind 

[Fe/H] = -0.116 ci - 1.624c? + 8.955 ci(fe-y) 

+ 42.008 (b-y)- 99.596 (b - yf + 64.245 (b - y)^ 

+ 8.928 ci mi + 17.275 nii - 48.106 w? 

-H 45.802m? -8.467, (3) 

which applies to stars with 0.43 < (b - y) < 0.63, 0.07 < nii < 
0.68 and 0.16 < ci < 0.49 with a standard deviation of 0.12 dex 
(the same cr is obtained considering instead Eq.|2]for equally red 
stars, but averaging with this latter form helps to reduce the ze- 
ropoint offset for cool stars) . With respect to the functional form 
used in Schuster & Nisse 3 ([1989), ours has the same standard 
deviation, but performs significantly better for [Fe/H] < -1 .0. 

The right hand panel of Fig. |4] shows predictions using 
"MARCS-standard" synthetic colours: models capture the main 
trends, especially at higher T^ff and different surface gravities, 
where the choice of various log g in the synthetic spectra is ap- 
proximated by the 10 and 90 percentiles of the ci distribution 
in the data (assuming lower values of ci to trace higher log^, 
which does not hold exactly towards the coolest Teff , because of 
contamination between dwarfs and subgiants). The main point 
from Fig. |4]is that our adopted functional form is a good rep- 
resentation of the data, even in poorly sampled regions of the 
plot and the trend at super-solar metallicities, where we do have 
calibration stars, is real. 

We applied Eq. ^ to all stars in the GCS, but for stars red- 
der than (b-y) > 0.43 we also used Eq. (O and then took the 
average of both estimates as our final value. The comparison be- 
tween the input spectroscopic metallicities and our photometri- 
cally derived values is shown in Fig.|5] Both equations provide a 
good representation of spectroscopic measurements and, within 
their accuracy, we do not introduce any obvious discontinuity. 
Our procedure gives a more homogeneous sample, avoiding the 
presence of breaks in different colour ranges, as was the case in 
the previous GCSII (see Fig. |6]l. Uncertainties in the observed 
Stromgren colours also bear on derived metallicities. On aver- 
age, the effect amounts to 0.04 - 0.05 dex in [M/H] and [Fe/H], 
estimated running a MonteCarlo simul ation with ob servational 
errors in (b - y), mi and ci as given in lOlsenI ([T98l . Errors in 
the derived metallicities tend to increase towards the blue- and 
red-most indices. 

A comparison with the homogenized spectroscopic sample 
of ' Arnadottir et al] (1201 Ol) confirms the quality of our calibra- 
tion with a median (mean) difference (ours minus Arnadottir) 
of 0.002 (0.007) dex and a scatter cr = 0.13 dex. Note that 
as discussed throughout the text, the overall scatter of our cal- 
ibration with respect to the spectroscopic sample is slightly be- 
low 0.1 dex, though this comparison folds the uncertainties that 
affect spectroscopic estimates. The test on open clusters (see 
Section 12.2.1b suggests that the intrinsic scatter in the metal- 
rich regime is actually somewhat lower. Usi ng the recently deter- 
mined uvby solar colours of Mel endez et al . (2010b), we obtain 
[Fe/H]Q = -0.006 dex, which agrees well with the zeropoint of 
our metallicity calibration. 

2.2.1 . Further test of the [Fe/H] scale 

We already checked our metallicity scale against the homog- 
enized stellar sample of lArnadottir et al.l (120101) and the solar 
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Fig. 6. A[Fe/H] (in the sense ours minus GCSII) as function of 
Teff (panel a) and [Fe/H] (panel b) for stars in the irfin sample 
satisfying the applicability range of our metallicity calibration. 
Notice the two breaks at 5500 K and 6500 K which correspond 
to the discontinuities introduced by the three calibrations used in 
GCSII for different (b - y). Contour levels and lines as in Fig.|2] 



colours of 'Melen dez et al ] (12010b'). Here we further test it by 
using open clusters and a moving group; finally we comment 
upon the limit of our calibration for intrinsically bright stars. 

Hyades and Coma are two nearby, virtually redd ening-free 
clusters often used to check the meta llicity scale (e.g lHavwoodI 
12006). though'Holmberg et al.' ( 2007) claim that GCS uvby pho- 
tometry of stars belonging to the Hyades cluster is not on the 
same scale as the rest of the catalogue, possibly because they 
were observed at higher air masses from Chile. In the case of 
Hyades, its controversial ci colour anomaly is also of concern 
(i.e. the systematic difference in the ci vs. (b - y) diagram be- 
tween the sequence of unevolved stars in Hyades and the cor- 
responding sequence for unevolved Coma and field stars w ith 
similar 6mi; see e.g.. lCrawfordlll97"5tlStromgren et al.lll982[ for 
more details). 

We took Strom gren photometry for the Hyades cluster from 
ICrawford & Perrvkil96 6 and compared the result of our metal- 
licity cahbration with the detailed spectroscopic study of Hyades 



' The extension of the original uvby system to cool a nd metal-poor 
stars is based on two main sets of standard stars, those of lBondl ( 119801) 
and Olsen ( 1993), respectively. The main discrepancy between the two 
concer ns the Ci inde x, stemming from dilferences in u band (see discus- 
sion in l01senlll995h . Our metallicity calibr ation uses GCS photom etry, 
which is built on the Olsen standards (see iNordstrom et"al]|2004l . and 
references therein). For testing in the metal-rich regime, as we do here. 



stars of iPaulson et al.l (l2003h (who adopt a T^g scale rather close 
to our IRFM scale). There are 10 single stars in common and 
we find a mean [Fe/H] = 0.09 + 0.02 dex (cr = 0.06 dex). 
This value is shghtly lower than the mean obtained using the 
same stars in lPaulson et al.l (l2003l) . which amounts to [Fe/H] = 
0.14 + 0.01 dex (cr = 0.04 dex) and nearly coincides with the 
mean value derived by ^ Paulson et all (120031) using a larger sam- 
ple of cluster members. We note that for this cluster a typical 
meta llicity around 0.1 dex is commonly cited in the literature 
(e.g.. lTavlor & Joneill2004HSchuler et all2006l) . 

The Hyades open cluster is known to be underabundant in 
helium for its metallicity (AY ~ 0.02 see e.g.. lVande nBerg et al.l 
2010. and references therein), a feature which would be tempt- 
ing t o associate to the ci colour anomaly dStromgren et al.l 
'1982"). However, synthetic colours show that variations of he- 
lium of this order affect ci to a negligible extent (Melende z et al.l 
2010b). Another possibility is that the anomaly is caused by 
variations in other elements. In fact the Hyades anomaly could 
simply be the [Fe/H] diff'erence between stars of similar 6mi 
as the following comparison with the Coma cluster suggests. 
For 17 stars in the Coma cluster we took the photometry of 
ICrawford & BarnesI (Il969l) and derived [Fe/H] = -0.08 ± 0.02 
(cr = 0.07 dex) using our calibration, which implies a metallic- 
ity diff'erence with respect to Hyades that excellently agrees with 
that spectroscopically measured by iBoesgaard & Friei (1 19901) 
and Friel &_Boesgaard ( 1992). 

Another open cluster originally observed by 
ICrawford & BarnesI (1 19701) is NGC752. This cluster 
also has relatively low r e ddenin g E(b - y) - 0.027 
(lAnthony-Twarog & TwarogI |2006|) . Using all dwarfs in 
ICrawford & Barnes (1970) (within the colour range of our cali- 
bration) gives |Fe/H| = -0.07 + 0.02 dex (cr = 0.1 1 dex), and a 
similar value (-0.05 ± 0.04 dex, cr = 0.09 dex) when restricting 
the same photometric measureme nts to the smaller -yet with 
cleaner membership- sample of lAnthonv-Twarog & TwarogI 
(12006 !). 

Finally, using observations of F-type stars in the Pleiades 
(which are less affected by activity stemming from the young 
age of this open cluster) from Cra wford & PeiTV (1976) and 
adopting E(B - V) — 0.04 (e.g., Ivan Leeuwenll2009l) . we de- 
rive [Fe/H] = 0.00 ± 0.02 dex (cr = 0.10). The diff^erence 
with respect to the Hyades again excellently agrees with that ob- 
tained from the spectroscopic comparison of Boesgaard & Frie| 
(1990), after co rrecting the Pleiades for known non-members 



( An et al.ll2007h . Our [Fe/H] also agrees well with recent spec- 
troscopic estimates based on a Tgif scale consistent with our 
own (Soderblomet al. 2009). For the last two clusters, we also 
checked that a typical uncertainty of E(b - y) - 0.01 affects 
[Fe/H] by ~ 0.01 dex. 

An additional check on the precision of our met allicity 
calibrations comes from the HR1614 moving group (Eggenj 
1978; Feltzing & Holmberg 2000). Chemical tagging via high- 
resolution spectroscopy of kinematically selected members al- 
lows us to clearly identify interlopers amongst the group mem- 
bers (iDe Silva et al.ll2007l) . Fig.|7]shows the differential Fe abun- 
dance A[ Fe/H ] for a number of candidate members in common 
between lOe Sil va et al. (2007) and GCS, using our metallicity 
calibration. The plot is relative to the mean metallicity of the 
sample, and thus largely independent on the underlying Tgff scale 
a dopted. The comparis on agrees remarkably well with figure 2 
in lDe Silva et alJ (l2007l) . clearly allowing us to identify spurious 



these differences are of no concern since the original set of observations 
defining the uvby system (and therefore adopted also by Olsen) is used. 
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Fig. 7. Differential Fe abundances of kinematically selected 
members of HR1614. Filled circles are metal-rich members, 
while asterisk (HIP 135 13), triangle (HIP6762) and diamond 
(HIP25840) are spurious members according to the chemi- 
cal tagging performed by iDe Silva et al.l (l2007b . Open circles 
are the analogous comparis on with the group members from 
iFeltzing & HolmbergldioOOl) . 



members of the moving group. We determine the group to have 
a mean [Fe/H] - 0.28 + 0.02 dex (cr - 0.07 dex), in good agree- 
ment w ith the spectroscopic value of 0.25 dex in lDe Silva et alj 
(12007'). 

Finally, we comment on the accuracy of our photomet- 
ric metallicities for intrinsically bright stars. The spectroscopic 
sample upon which our calibration is built extends to magnitudes 
only slightly brighter than My-^ ~ 2, which are typical for F 
dwarfs; however, the GCS contains some hundreds of stars more 
luminous than this (also compare with Fig. [T2l i. These stars are 
close to the instability strip and are therefore possibly contami- 
nated by 6 Scuti pulsators and/or chemically peculiar A/F stars 
(the latter often bei ng overabundant in Fe and possibly w ith pe- 
cuHar colours, e.g. jGebran et al.ll2008HNetopil et al.ll2008 ). Our 
calibrations include a dependence on the ci index (a good sur- 
face gravity indicator for hot stars), so in principle, we can ex- 
pect them to work for decreasing log g. Fig.[8]shows a clear trend 
for the brightest stars in the GCS, which tend to be more metal- 
rich than the remaining part of the sample (left panel). Because 
they are preferentially metal-rich, an age determination based 
on isochrones also biases them to even younger values. These 
bright stars also stand out in the study of the metallicity distri- 
bution function (Section|4]i and in kinematic (Fig.fTTTi. Some up- 
trend in this figure (starting around My^ ~ 3) is expected from 
colour/spectral-type cuts in the original GCS sample selection, 
because at a given colour metal-poor dwarfs are fainter (cf. also 
Fig.O. However, at the bright end (in particular from Mvj ~ 2) 
the calibration seems to deviate too strongly. Interestingly, those 
stars are preferentially the most distant ones and thus have in- 
creasing reddening uncertainties as well as the largest paral- 
lax errors, which could misleadingly place intrinsically bright 
stars at fainter absolute magnitudes. Also notice that by sam- 
pling larger distances, where the GCS is not complete anymore, 
intrinsically luminous stars are preferentially found around the 
peak of the metallicity distribution function, which could partly 



account for their rather high [Fe/H]. In addition, preferentially 
higher metallicities could also stem from these objects originat- 
ing from the inner disc (cf. Section |6), a feature of which we do 
not find any significant indication, though. Even if not conclu- 
sive, this seems to suggest that removing stars with an absolute 
magnitude Mvj- brighter than 2 would be a safe choice when us- 
ing GCS stars for deriving local constraint on Galactic chemical 
evolution. 

2.3. The mild sensitivity of Stromgren pliotometry to fhe 
a-eiements 

For all 1498 calibrating stars presented in Section 122] we have 
[Fe/H] and [a/Fe] from high-resolution spectroscopy. From this 
the o verall m etal-to-hydrogen ratio [M/H] can be computed 
("e.g. jYi et al.ll200 1). Interestingly, when fitting functional forms 
of the kind of Eq. and Q to [M/H], the scatter of the re- 
sulting calibrations decre ases to 0.07 and .10 dex respectively, 
thus suggesting (cf. e.g.. lYong et al.ll2008h that Stromgren in- 
dices carry information on the overall metal content (apart from 
a few of the most metal-poor stars in the sample, which indi- 
cates diminishing sensitivity to [M/H] because of the intrinsi- 
cally fewer lines). 

It would be possible to apply our calibrations as a func- 
tion of [Fe/H] and [M/H] and from those derive an estimate of 
[ff/Fe]. In practice though, there is some degree of correlation in 
the results since the same functional form and indices are used 
over 2 dex in metallicity to estimate typical alpha-enhancements 
within ~ 0.5 dex. We experimented with different combinations 
of Strorngren colours and found ai - {v - y) - (b - u) to be 
sensitivqjto [a/Fe] at a given T^g. Fig. |9] shows [a/Fe] versus 
our index ai as well as the comparison with synthetic colours 
at a few (b - y) values for the sets of alpha-enhanced and -poor 
models available through the MARCS library. 

In Fig.|9]a dependence on [Fe/H] is certainly built in given 
that stars with lower [Fe/H] have preferentially higher levels 
of alpha-enhancements. Nevertheless, the comparison with syn- 
thetic colours shows that the trend is real at fixed metallicities. 
As we akeady pointed out, Stromgren synthetic colours are not 
immune to deficiencies, and combinations of the adopted filters 
are -by construction- also sensitive to metallicity and surface 
gravity. Changing the latter parameter shifts synthetic colours to 
the right or left with respect to the position shown in Fig. |9] 
which refers to logg = 4.5. However, we checked that the 
shape of the slopes remains unaffected by the exact value of 
logg. Limitations in synthetic colours as well as surface grav- 
ity dependence could explain why the bulk of calibrating stars 
is fitted by models having [Fe/H] - -0.5 rather than a higher 
metallicity, which is more representative of the sample (cf. also 
with Fig.|4]i. According to the models, in Fig. |9] the sensitivity 
to [a/Fe] is more pronounced (i.e. it has a shallower slope) at 
cooler effective temperatures, which are therefore likely to be 
better recovered. Determining [a/Fe] becomes increasingly dif- 
ficult for the hottest and most metal-poor stars, as expected be- 
caus e both atomic and m olecular lines get weaker in this regime 
(e.g. ICoelho et al.ll2 005). Yet, even at the bluest colours the data 
seem to show a clearer trend with [(r/Fe] than models. Aware 
of these warnings, the mild correlation of the oi index with 
alpha-enhancement seems to work for drawing meaningful con- 



Also other indices have been found to show some dependence on 
[a/Fe] such as e.g., mi - (b - y). From our investigation it seems that 
Stromgren filters such as b and y are barely affected by [a/Fe], while ii 
and V are more affected. 
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elusions wfien one fias a statistically large sample of stars (see 
also FigfTTli. 

For eacfi (b-y) in Fig.|9]we constructed a fiducial using stars 
of similar Tetj and derived a value of [a/Fe] according to tfieir 
(v - y) - (b - u) witfi respect to tfiat of tfie corresponding fidu- 
cial. Despite models sfiow a spread witfi metallicity in Fig. |9] 
we did not include any dependence on [Fe/H] in building tfie 
fiducial to avoid any risk of introducing a spurious trend of in- 
creasing alpfia witfi decreasing [Fe/H]. Tfie comparison between 
tfie spectroscopic measurements and our pfiotometric estimates 
is sfiown in Fig. [10] Tfie overall agreement is indeed good (for- 



mally cr = 0.09 dex), tfiougfi tfiere are a few caveats: [a/Fe] for 
stars witfi [Fe/H] < -1 is not well recovered (tfie calibration sat- 
urates, filled squares), as expected from our previous discussion 
on metal-poor stars. Also, for tfiin-disc stars [a/Fe] tends to be 
sligfitly underestimated/overestimated at fiigfier/lower metallici- 
ties. Tfiis reflects tfie sfiape of tfie fiducials used to derive [a/Fe] 
from Fig. |9] In fact. Fig. [TT] sfiows tfiat our [a/Fe] calibration 
does not allow us to recover any gap between tfiin and tfiick disc 
stars. Tfie sfiape of tfie overall narrow trend is tfius driven from 
tfie fiducial, yet witfiin tfiis trend a distinction between alpfia- 
ricfi and -poor stars is possible, tfiougfi in a statistical sense only. 
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Fig. 10. Upper panel: spectroscopic versus Stromgren [a/Fe] es- 
timates for our 1498 calibrating stars. Crosses (circles) are stars 
having probability > 90 percent of being thin (thick) disc based 
on their kinematic (using U, V, W velocities from the GCS). 
Downward triangles are stars with lower probability or for which 
kinematic information was not available. Filled squares are stars 
having [Fe/H] < -1 (independently of their kinematic thin/thick 
membership, if available). Lower panel: same symbols as above, 
showing the difference spectroscopic minus ours. 



This is shown by selecting calibration stars on the right (left) 
of the dotted (dashed) line in Fig. [TT] with the same distinction 
right vs. left still being preserved when we use our [Fe /H] and 
[o-ZFe] calibrations, which are represented by upward vs. down- 
ward triangles. Notice that using our [a/Fe] the dispersion of 
the photometric [M/H] with respect to the spectroscopic mea- 
surements is 0.08 dex, compared to 0.10 dex (previous section) 
when using [Fe/H] only. 

In conclusion, the sensitivity of our approach to the alpha 
elements is real, but mild and works only for [Fe/H] > -1 
or slightly higher values. Also, a statistical distinction between 
alpha-rich and -poor stars is possible, but only within the func- 
tional form of our calibration so that other finer structures could 
still be missing. Thus, the values of alpha-enhancements we de- 
rive are not exact measurements of [(r/Fe], but rather a proxy 
of them for stars of similar [Fe/H]. For this reason we will re- 
fer to our estimate as orFe instead of [a/Fe] throughout the pa- 
per. Nevertheless, as we show below, when one deals with sev- 
eral thousands of stars, as is the case in the GCS, our cFe can 
give important insight into the formation and evolution of the 
Galactic disc(s). 

3. New age and mass determinations 

Revising metallicities and effective temperatures also affects age 
and mass estimates for the stars. Fig.[T2]shows comparisons be- 
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Fig. 11. [Fe/H] versus [a/Fe] for the 1498 stars in the spectro- 
scopic sample (filled circles). Dotted and dashed lines are arbi- 
trarily used to select stars on the right and left locus of the plot 
(only for [Fe/H] > -1). Downward (upward) triangles represent 
stars selected according to this criterion, but using our [Fe/H] 
and [a/Fe] Stromgren calibrations. 



tween isochrones and stars with metallicities close to those of 
the plotted isochrones. Compared to previous studies, our im- 
proved effective temperatures are hotter, and the large system- 
atic discrepancies betw een theoretic al isochrones and observed 
data that plagued e.g., Pont & Eyed ([2004.) almost entirely dis- 
appear. As can be seen, only at the lowest metallicities and lu- 
minosities the theoretical main sequence has the tendency to fall 
beneath the stars, i.e., the isochrones are too hot. However, the 
discrepancy is considerably reduced from earlier GCS analyses 
where shifts in the effective t emperature of the isoch rones have 
to b e introduced below so lar (iNordstrom et al.ll2004b or even at 
all dHolm berg et al.l l2007h metallicitiefl Our new temperature 
and metallicity scales prove to agree very well with those of the- 
oretical isochrones, at least for metallicities higher than about 
[M/H] > -0.5, which includes the vast majority of stars in our 
sample (cf. Fig. [TSl l. Indeed, for the sake of Fig. [12] a metallic- 
ity bias stemming from the wings of the metallicity distribution 
function will also play a role, as we discuss in greater detail in 
the appendix, as well as the monotonic decrease (increase) of 
stars with metallicity in the metal-rich (-poor) tail of the metal- 
licity distribution function by which the average metallicity in 
a given interval can be lower (higher) than the middle value of 
the interval. This can be clearly seen in the net bias of hotter and 
more metal-poor (cooler and more metal-rich) stars in the top 
left (bottom right) panel of Fig. [T2] where stars are selected in 
symmetric intervals around the metallicity of the isochronej^ 



' Note that in GCSII also the solar isochrone, which is constructed to 
fit the Sun, has to be cooled by 0.005 dex in log Teff, corresponding to 
approximately 70 K, in agreement with the ofi'set in Fig. [2] 

At the same time, a stronger disagreement for metal-poor low main 
sequence stars -for which their position on the HR diagram is substan- 
tially age independent- was noticed by Casagrande et al. (2007) and it 
could have potential implications for studies of multiple stellar popula- 
tions (Portinari et al. 2010). The same Teff scale adopted here compares 
well with isochrones for nearby, evolved subdwarfs, suggesting that this 
disagreement s eems now reduced at lea st for log Teff > 3.7 (compare 
with fig. 10 in IVandenBerg et ani201d) . though further investigations 
are encouraged. 
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Fig. 12. BASTI isochrones for different ages at a given metallicity (continuous lines) compared to stars of similar [M/H], where the 
difference ±0.1 or +0.2 dex is coded by colours. Larger symbols are for stars with higher parallax accuracy as labelled in the top 
left panel. Only stars in the irftn sample are shown. 



As pointed out by iPont & Eyed (l2004l) . naive fits to 
isochrones lead to severe biases, e.g. what they name a termi- 
nal age bias. This happens because some places on isochrones 
are more densely populated than others because of the map- 
ping from mass to colours/luminosity owing to the initial mass 
function and to the time scales involved in stellar evolution. Just 
looking for the closest match ignores these facts and might er- 
roneously place too many stars into sparsely populated regions. 
Biases of this kind can be accoun ted f or by taking a Bayesian 
approa ch as in iPont & Ever! (1200 4^) and ^J 0rgensen & Li nde grenI 
(l200l . who did a Bayesian age determinatio n on the old GCSI- 
II. A detailed discussion can also be found in Burnett & BinnevI 
(1201 0). In our sample the errors vary significantly between stars, 
but they depend only weakly on the derived stellar parameters, 
so that we can neglect this influence on the age distribution. We 
only used log(reff), absolute Johnson V magnitude and metal- 
licity information to estimate the ages and masses of stars. In 
principle more information could be in the colours, but essen- 
tially this is already exploited by the colour-dependent calibra- 
tions. Moreover, a direct use of colour information would im- 
ply relying directly on synthetic colours (with th e uvby bands 
being more troublesome than others, see e.g., iOnehagetalJ 
120091: IArnad6ttiret"ani2010l: iMelendez et 311120 lObh . which we 
wanted to avoid. Further, especially Stromgren colours depend 



on [M/H], thus requiring an even denser grid of metallicities 
than the one we use. 

The isochrones provide us with a natural grid for calculating 
the probability distribution function for the parameters of a given 
star Every isochrone point has to be weighed by the volume of 
parameter space it has to cover and by the a priori assumptions. 
To avoid any factor that could contribute to the particular age 
distribution that we find in the sample, we assume a flat age prior 
for 0-14 Gyr, i.e. a constant density J?[(t) of stars over age. 

[1 for < T < 14Gyr 

[O else 

As different positions of stars in the Hertzsprung-Russell dia- 
gram imply different underlying selection functions that bias the 
intrinsic age distribution at this place, this approach also avoids 
making further assumptions that could potentially weaken the 
interpretation o f the results. The mass prior is a Salpeter IMF 
(ISalDeteiill955h and we do not set any dependence on age. In 
the mass interval of interest here (cf. Fig. [14), a Salpeter IMF 
is indeed still appropriate, whereas considerably larger uncer- 
tainti es exist regarding the lower and higher mass range (e.g., 
Bas tian et a D|2010). We further tested the extreme case of a flat 
IMF, and even this unrealistic assumption has negligible impact 
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Fig. 13. Normalized age probability distribution for all stars in 
the sample having both BASTI (continuous thick) and Padova 
(continuous thin) ages and only for stars in the irfm sample (dot- 
ted thick and thin lines) having good ages and Mvj > 2. 



on our results (see the appendix, where we provide details on 
the Bayesian scheme adopted for dealing with the observational 
errors in T^ff, metallicities and absolute magnitudes). 

In order to study differences between different isochrones, 
we used grids of the BASTO (Pietrinfemi et al. 2004, 200^ 
^9) and PadovS isochrones (BertelH et al. 2008, 2009). The 
Padova grid has a logarithmic age spacing of 0.01 dex, i.e. it 
rises from 23 Myr at r = 1 Gyr to 230 Myr at t = 10 Gyr. 
We queried a total of 56 metallicities from the database, which 
are created by interpolating among the nine available metal- 
licities, ranging from Z - 0.0001 to Z - 0.07. The solar 
isochrone has Yq - 0.26 and Zq - 0.017 (iGrevesse & Sauvall 
Il998l) . The helium-to-metallicity enrichment ratio was chosen 
to be AF/AZ = 2.1, which is consistent with the value inferred 
from the study of metal -rich local K dwarfs (Casagrande et al. 
I2OO7I) . At the lowest metallicities this falls somewhat short of 
Y - 0.23 (the lower helium abundance in the database, lower 
than the curren t preferred estimate from WMAP +BBN, see e.g., 
IStei gman 2010). and for those objects we kept this value of Y. 
For the B ASTI iso chron es (Fq = 0.2734 and Zq = 0.0198 from 
IGrevesse & Noelsl 1 1 993h we used a denser grid than the pub- 
lished one. This grid was specially calculated for this purpose 
to include 20 metallicities at AF/AZ = 1.45 (leading to a pri- 
mordial helium abundance in agreement with the cosmological 
estimate) in the range Z = 0.0001 to 0.04, with a time spac- 
ing of 100 Myr maximum, making this grid denser (sparser) at 
high (lower) ages compared to the logarithmic age spacing in the 
Padova isochrones. Both sets of isochrones assume solar-scaled 
abundances (i.e. constant ratio of the single metals with respect 
to the Sun), which are appropriate because it has been shown 
that for the range of metallicities covered by the present study, 
isochrones for enhanced a abundances can be reproduced re- 
mark ably well by those for solar scaled mixtures if Z is the same 
(e^.. 'Chieffi et al."1991': 'Chabo ver et al.lll992l; ISaiaris & WeissI 
11998; VandenBerg et al, 2000). We also checked the difference 
when using [Fe/H] rather than [M/H] in determining ages: the 
overall difference is fairly small, with a scatter of about 0.5 Gyr. 



' ' http://www.oa-teramo.inaf.it/BASTI 

http://stev.oapd.inaf.it/YZVAR/cgi-bin/foi-m 



Fig. [13] shows the age probability distributions for all stars 
in the GCS with ages determined from both BASTI and Padova 
isochrones and also for all stars with good ages. Throughout the 
paper, ages are defined to be good if cr < 1 Gyr or the rela- 
tive uncertainty is better than 25 percent (see also the appendix). 
While these criteria are arbitrary, they balance a reasonable de- 
termination of absolute ages for young objects with a reasonable 
relative determination for older ones. The distribution strongly 
peaks around 2 Gyr, which is caused by th e selection effects on 
the sample (see also lNordstrom et ani2004l) . 

The GCS is in fact limited near the plane of the disc, while 
older stars usually have a considerably more extended verti- 
cal distribution, which brings their orbits high above the plane 
and thereby lowers their presence in this survey. In addition, 
the magnitude limits of the catalogue give a larger volume to 
bright, young stars, and the exclusion of giant and very blue 
stars from the sample leads again to a net bias against very 
young and especially against old objects. An estimate of the 
age of the disc can thus not be done directly using the age 
of the stars in the present sample, but requires modelling the 
star-formation history of the solar neighbour hood, returning a 
considerably older disc (> 10.5 G yr, see e.g., lAumer & BinnevI 
120091: ISchonrich & BinnevI l2009al) . Also notice that because of 
the young ages, our sample is fairly immune to atomic diffusion, 
possibly apart for a few of the oldest stars. 

Fig. [T4| clearly summarizes all main issues in dating stars. 
Ages are most readily determined for stars in the upper en- 
velope of Fig. [141 which roughly maps the turn-off region. At 
low masses, apart from the most metal-poor subdwarfs, reliable 
ages are difficult to derive (grey dots) because the majority of 
these stars are still on the main sequence due to their long life- 
times. In addition, somewhat below 1 Mq the GCS starts losing 
completeness, being mostly limited to FG dwarfs. The youngest 
stars cover a short mass range (cf. blue points in the middle 
panel of Fig.[T6li: more massive young stars are in fact brighter 
and hotter than sample selection limits, apart from a handful of 
bright objects (squares). The reliability of our metallicity cali- 
bration for those stars was already discussed in Section 12.2.1! 
The depletion of stars longward of the kink at ~ 1.5 Mq pre- 
cisely suggests that at masses higher than this value the sample 
is partly incomplete (which roughly corres ponds to My^ ~ 2, 
using the mass-luminosity rel ation of e.g. iHenrv & McCarthvl 
ll993UFang & Yan-Nin"gll2010l) . On the contrary, no obvious bi- 
ases seems to be present in the range 1.1 < < 1 .5. 

4. The metallicity distribution function 

Given its complete nature (see Section|2]i, the GCS is well suited 
for the study of the met allicity distribution function (MDF) in 
the so lar neighbourhood (iNordstrom et al.l2004i iHolmberg et al.l 
^007"). However, this does not mean that the MDF shown here 
can be directly compared to theoretical expectations. For the 
same reasons presented above when discussing the age distribu- 
tion, sample selection effects enter the results. For a quantitative 
comparison those selecti on effects have to be taken in to account 
in theoretical models (cf. lSchonrich & Binnev"2009a!). 

We already argued in Section [2] that dividing the original 
sample into two groups does not introduce any bias. Stars with 
the best photometry show lower dispersion, but the average prop- 
erties are robust and are the same for both the irfm and clbr sam- 
ples. This is shown in Fig. [15] for [Fe/H] and [M/H], with the 
relevant statistical parameters given in Table [T] A Kolmogorov- 
Smimov test between the irfm and clbr samples for [Fe/H] and 
[M/H] tells that the probability of both samples being drawn 
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Fig. 14. Ages versus masses for stars belonging to the irfm sam- 
ple. Colours are for stars with well determined ages, going from 
metal-poor (blue) to -rich (red), while grey dots are for the re- 
maining stars. Squares are stars brighter than My^ - 2. Inner 
panel: same as outer panel, but with a metallicity colour coding 
also for stars with less reliable ages. 
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Fig. 15. MDF of the solar neighbourhood in terms of [Fe/H] (up- 
per panel) and [M/H] (lower panel). Continuous line refers to 
stars belonging to the irfin sample (5976 stars within the colour 
ranges of the metallicity calibration), dashed line when consid- 
ering only stars fainter than My,, - 2, dotted line to all clbr 
stars (8470 within the colour ranges of the metallicity calibra- 
tion) and dot-dashed when applying the same luminosity cut as 
above. Poisson error bars are shown for a representative case in 
both panels. 

from the same distribution is below 1 percent, i.e. not significant. 
The reason for this lies in the broader wings of the clbr sample, 
partly because the lower quality of the latter sample could be 
responsible for less reliably determined metallicities that over- 
populate the wings, and/or older ages (see below). When restrict- 
ing the selection to -0.5 < [Fe/H] < 0.5, the irfm and clbr sam- 
ples are in fact drawn from the same distribution to a level better 
than 5 percent, under the null hypothesis that the two distribu- 
tion are drawn from the same parent population. Identical con- 
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Fig. 16. Top panel: MDF for stars belonging to the irfin sample 
divided into different age intervals. Stars having age < 1 Gyr are 
shown with a continuous line, 1 < age < 5 Gyr with a dashed 
line and age > 5 Gyr with a dot-dashed line. Shaded areas iden- 
tify the subgroup of stars in the same age intervals as above, but 
with absolute magnitudes (< 2); no such bright stars are present 
in the old sample. Only stars with well determined ages (see 
SectionO are used. Bars indicate Poisson errors. Middle panel: 
[Fe/H] versus stellar mass. Colours have the same meaning as 
in the top panel, with grey dots now referring to the remaining 
stars having more uncertain ages. Filled squares identify stars 
with bright absolute magnitudes (< 2). Lower panel: same sym- 
bols and colours as in the middle panel, but showing the age- 
metallicity relation. Shown for comparison (asterisks) are the 
ages and metallicities of the halo Globular Clusters studied in 
iVandenBerg et al. 201 (in the latter case, a different zeropoint 
on the age scale is possible, also depending on the input physics 
adopted in the stellar models employed). 

elusions to the Kolmogorov-Smirnov statistic are reached using 
instead the Wilcoxon Rank-Sum test for comparison. We find 
that the MDF for young and old stars look considerably different 
(see below). We note that because the clbr sample contains a few 
more cooler stars than the irfin sample (see Section 12.1.21 ). the 
cooler stars being preferentially older and thus with a broader 
MDF (see below), this could also be partly responsible for the 
different broadening of the wings. 
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Slicing the MDF into different age intervals shows an inter- 
esting feature: young stars have a considerably narrower distri- 
bution than old stars, though the peak always remains around 
the solar value (Fig. \T6[ . Notice that because of the selec- 
tion effects on the sample age, an uneven slicing -denser at 
young ages- is more appropriate (cf. Fig. [13). While the MDF 
has been historically used to constrain the gas infall rate (e.^. 
Lynden-Bell 19 75; Jinsley 1980; Matteucci & Francois 198S 



Table 1. Metallicity distribution function 



Chiappini et al]ll997l) . the increasing broadening with age sug- 



gests that old stars are also a relevant ingredient in describ- 
ing the wings of the MDF. A n atural explanation is prov ided 
by the radial migration of stars (Sellwoo d & Binnevll2002l) . In 
this picture the solar neighbourhood is not only assembled from 
local stars, following a local age metallicity relation, but also 
from stars originating from the inner (more metal-rich) and 
outer (more metal-poor) Galactic disc that have migrated to 
the presen t position on different timescales (Roska r et al.li2008i: 
ISchonrich & Binnev. 2009a). Because of the higher density of 
stars in the inner disc, migration would favour metal-rich stars, 
which could compensate the metal-poor tail typical of local 
chemical evolution, which would explain the rather symmetric 
shape of the MDF we derived. A more quantitative explanation, 
however, requires modelling of the chemical evolution. 

The presence of a metal-ric h tail in Fig . [T6]could be a sig- 
nature of the Galactic bar (e.g.. ' Grenonl[T999ft : such a detection 
is however very difficult to claim even with the current sample. 
Indeed, we only detect a conspicuous young metal-rich popula- 
tion at the brightest magnitudes, where the accuracy of the metal- 
licity calibration could be lower (see the discussion in Section 
I2.2I 1. The presence of a bar would rather imply the existence of 
an old metal-rich populatio n, which we do no t detect (but see 
iMinchev & Famae v 2010; Minchev et al.ll20l"ll for a recent dis- 
cussion on the effect of the bar). Although we do not have ac- 
cess to the sample selection performed in the original assembly 
of the GCS, we regard the presence of a bias against old metal- 
rich stars as unlikely, and we refer to iNordstrom et al.l (1200 4l) 
for more details on the completeness of the sample. We also in- 
vestigated whether the metal-rich stars display any conspicuous 
feature in the UV velocity plane and did not find any. Notice 
though that the fraction of these young metal-rich stars in the 
total sample is fairly small and they do not bear considerably on 
the overall MDF of Fig. [B] 

Apart from the aforementioned bright stars, the metal-rich 
wing of the MDF is not an artefact caused by the sampl e se- 
lection on colours (contrary, e.g., to iKotoneva et al.ll2002h . be- 
cause high-metallicity stars are present throughout the entire 
mass range (middle panel in Fig. [T&t . Also, on the metal-poor 
side there is a clear contribution of (nearly) unevolved subd- 
warfs -for which a determination of ages is more uncertain- 
with a trend in mass mirroring that already observed in lumi- 
nosity (cf. Fig. [8] and[T4ll. 

The peak of the MDF is only slightly subsolar (me- 
dian rFe/H1 ~ -0.0 5 rM/H1 0.01), in a greement with 

Havwood I (120011). 'Tavlor & Croxall (l200l . Luck & Heited 
and Fuhrmann (2008), but in contrast with other 
studies, which rather favour a peak in the range -0.2 to 



(12001 " 



-0.1 de x (e.g., [W yse & Gilmore 1995; Rocha-Pinto & Macidj 
1 9961: lA llende Prieto et all |2004 INordstrom et al.l |206¥ 
Holmberg et al.. ,2007i) . In most cases the reason for this 
difference stems from the T^ff scale we use, which supports 
spectroscopic studies that adopt similar effective temperatures 
and results in higher metallicities. As a si de remark, we note 
that the MDF de t ermiri e d from M d\yarfs (iBonfils et all 120051: 
ICasagrande et all 120081: ICasagrandel l2008l) is hkely to peak 
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Notes. Statistical peak values of the MDFs of Fig.[T5]using stars in 
the irfm and clhr sample. Notice that the MDF is influenced by a 
low-metallicity tail. A Gaussian is not its best description. Median 
and FWHM provide different -and formally better- estimates. 
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Fig. 17. Velocity dispersion .v = ^(rfj + • ^y^, 
age. Age probability distribution functions derived from Padova 
and BASTI isochrones are used. Dotted black lines are 1 o" errors 
for the black line. In all cases stars with Mvj- < 2 are excluded, 
their effect being responsible for the bump (green dotted line) 
around 1 Gyr. 



aro und solar metall ic ity if the recent spectroscopic findings 
of ' Johnson & AppsI ( |2009|) are confirmed and photometric 
determinations for those stars are recalibrated accordingly. 

The peak at nearly solar metallicity of the local MDF at 
all ages also has implications for understanding secular pro- 
cesses associated with disc evolution, by i nvestigating whither 
and whence the Sun is moving (e.g., IWielen et al.1 119961; 
Bland- Hawthorn et a and it is also an im portant test of 

the overall solar metallicity (lAsplund et al.ll2009l) . 



5. The age-dispersion relation 

Figure [17] shows the velocity dispersion s for stars in the irfin 
sample as a function of stellar age. Ages are determined using 
the BASTI isochrones, apart from one case where the result of 
using Padova isochrones is shown for comparison. The differ- 
ence between requiring well determined ages (according to the 
definition of Section[3] black line) or not (cyan line) suggests that 
the signature of a continuous rise becomes even more promi- 
nent, confirming earlier stu dies of the GCS (Nor dstrom et all 
l2004t iHolmberg et al.l 120071) . This has to be expected; because 
of the pronounced overdensity of stars around ages of 2 Gyr 
(cf. Fig. [T3T l. excluding unreliable ages gives less contamina- 
tion to the rarer very young and especially to the older stars. 
However, because velocity dispersion roughly incr eases with age 
to the power 1/3 (ISpitzer & Schwarzschildl[l953l) . and because 
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Fig. 18. Panel a) and b): metallicity as function of orbital and 
guiding centre radius, respectively. The probability of a star to 
belong to the thin disc is represented by colour Stars with halo 
membership higher than 50% are plotted in black. Filled dia- 
monds and squares identify stars having Vlsr < -40km/s and 
> 20km/s, respectively (cf. with Fig. |20]|. Panel c): cumulative 
metallicity gradient ([Fe/H] and [M/H]) when including stars 
of increasingly older ages. Error bars are shown in one repre- 
sentative case. The gradient is computed using the mean orbital 
radius of stars as baseline. A kinematic cut to exclude halo stars 
(as described in the text) is also adopted for comparison. Panel 
d): same panel c), but using the guiding centre radius as baseline. 
Panel e): metallicity gradient centred at different ages, weighting 
all other stars with a Gaussian of width 1 .5 Gyr and using the 
mean orbital radius as baseline. Panel f): same as panel e), but 
using the guiding centre radius. Only stars with well determined 
ages (Section[3l) are used in all instances. 



of the ~ Gyr uncertainty in ages, it is actually difficult to distin- 
guish between a plateau and a real increase. 

When no metallicity nor kinematic cut is applied, a strong 
rise appears at the oldest ages (blue line). This feature is likely 
caused by contamination of moderately metal-poor stars that 
might belong to the Galactic halo. This disappears when us- 
ing a very conservative cut at [Fe/H] > -0.5 dex or a milder 
one at [Fe/H] > -0.8 dex but only considering stars with V > 
-150km/s. These cuts exclude some tens of stars, consistent 
with expectations from local disc-to-halo normalization which, 
despite lar ge uncertainties, is in the range of a few hundred s- 
to-one (e.g. 'Momson"l99l iGould et al.l[l998t lJuric et al.l l2608l) 
Because the isochrones might fail to exactly match metal-poor 
stars (Section O, the derived age distribution of low metallic- 
ity stars can be biased to older ages. Difficulties in understand- 



ing selecti on criteria are hkely to be responsible for the different 
findings of Ouillen & Garnett (2000) who -essentially using the 
sample of 189 stars studied in Edvardsson et al. ( 1993)- claimed 
the presence of a plateau in the dispersion over all intermediate 
ages followed by a quick rise at about ~ 10 Gyr 



6. Disc 

6. 1 . Metallicity gradient 

Abundance gradients across the Galactic stellar disc provide fun- 
damental constraints on the chemical evolution of this compo- 
nent of the Milky Way, and on the phy sical assumptions adopted 
in chemical evolutio n models (e.g., iPortinari & Chiosil 119991 : 
IChiappini et al.ll200lT) . Despite its local nature, the large number 
of stars in the GCS would suggest that it is possible to use it for 
estimating the radial m etallicity gradient in the Galaxy (cf. e.g., 
iNordstrom et al.ll2004 . 

While Galactic radial positions (i^cai) are snapshots of stars 
at the present time, covering a very limited range in distances (at 
most 0.3 - 0.4 kpc for the GCS), their mean orbital radii (left 
panels in Fig.fTSll allow us to probe larger distances (up to a few 
kpc) and thus are better suited for deriving the rnetalli city gra- 
dient (e.g. lNordstrom et al .120041 iHolmberg et al.l2007b . Orbital 
radii depend on the adopted Galactic potential; Rcai and V ve- 
locities offer an alternative and model-independent approach via 



the guiding centre radius / 



«G.,i(V+232) 
220 



under the assumption of 



a constant circular rotation speed of 220km/s (right panels in 
Fig.m. 

Using only stars that belong to the irfin sample and are within 
the calibration range, the exact value of the gradient still depends 
on whether or not a cut at the lowest metallicities is imposed 
to exclude contamination from halo stars. In Fig. [18] we show 
the case of applying neither kinematic nor metallicity cuts, as 
well as a kinematic selection to retain only stars with probabil- 
ity high er than 90 percent of belonging to the thin or thick disc 
(see e.g. lRamfrez et al ]l2007i for more details on this kind of se- 
lection procedure). Cutting the sample to exclude metallicities 
lower than about -0.8 dex has a similar effect as the kinematic 
selection. Metal poor (halo) stars having small orbital (and guid- 
ing centre) radii are in fact responsible for the strong positive 
rise in the metallicity gradient. 

It appears obvious that taking all stars at their face values 
does not provide a meaningful measure of the gradient in the 
disc. Indeed its value depends on the adopted kinematic or metal- 
licity cuts, the age interval considered, and also whether orbital 
or guiding centre radii are used for the computation (Fig. [T8j. 
Difficulties in estimating e.g. the interdependence between age 
and kinematic cuts (as stars with increasing asymmetric drift 
are preferentially older) as well as the increasing scatter in the 
age-metallicity relation and in the age-dispersion relation fur- 
ther complicate the picture. 

Fig. [TS] (middle panels) suggests the presence of a moderate 
negative radial gradient, consistent with studies using other indi- 
cators at various Galactocentric distances such as Cepheids, HII 
regions, B stars, open clusters and planetary nebulae (see e.g., 
Maciel & Costa 2010, and references therein). When restrict- 
ing the analysis to different age intervals (lower panels) there 
is an indication of a flattening and even a reversal of the gradient 
with increasing age, but we stress once more that the adopted 
kinematic or metallicity cuts affect the results. Such a signa- 
ture comes from older thick (as well as halo) stars (Spagna et aTl 
I2OIO), while the GCS is mostly limited to the younger objects 
that are situated in the thin rather than the thick disc. 
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As already mentioned in Section |4] the different behaviour 
for younger and older stars can be understood in terms of ra- 
dial migration, where the increasing age that is responsible for a 
broadening of the MDF could also soften the gradient, but more 
data and extended analyses are needed to explore this scenario. 

6.2. Thin, thick or stirred? 

Observations of external edge-on galaxies sh ow the presence 
of both a th in and thick disc component ('BursteinI 1979t 
iDalcanton & Bernstein 2002; Yoachim & Dalcanton 2006). The 
Milky Way seems to have a two-component disc as well, 
which was first propose d to fit the vertical density profile de - 
rived from star counts (lYoshiil Il982t iGilmore & ReidI Il983l) . 
Disentangling the nature and origin of these components is 
therefore highly relevant for understanding galaxy formation. 
While models in which thick and thin discs form sequentially 
via a rapid or dissipative collapse of proto galactic clouds be- 
came disfavoured during the past years (e.g.. lMaiewskil[T993l) . it 
is not yet clear how the stellar disc can form a thick component 
with time, if this is caused by to internal ( scattering, dynamical 
interaction or radial mixing, e.g., Scho nrich & Binnevll2009bl 
iLoebman et al.l 120101) or external (satellite accretion, mergers 
of gas-rich systems, minor mergers, e.g., IXba di et al. 20031 
Brook et al. 2007'; IVillalobos & Helmill2008t IScannapieco et^aP 
2009) mechanisms. 

Though limited to the solar neighbourhood, the GCS can 
provide important insights into this puzzle, because it is essen- 
tially free from kinematic selections. Our metallicities and cFe 
(Section [2.3b provide for the first time a way to investigate this 
with a more complete sample. 

Figure [19] shows all stars with a reliably determined [Fe/H] 
and Q-Fe and for which the U, V, W velocities are known, so that 
the same kinematic probabilistic selection scheme to the thin 
or thick disc adopte d in the previous section can be applied 
(iRamfrez et al ] l2007h . The small scatter and overall shape of 
the plot simply reflects the fiducial used to derive aFe, which 
squeezes up most of the metal-poor stars and also prevents us 
from seeing any gap between the thin and the thick discs (see 
discussion in Section IZJt . Despite these limitations, a qualita- 
tive picture can be drawn. Stars kinematically attributed to the 
thick disc populate the upper envelope of the Fig. [T9]for subso- 
lar metallicities, while merging into the thin disc around solar 
[Fe/H], in agreement with similar findings obtained by studies 
based on high-reso lution spectroscopy (e.g., iReddv et ani2006t 
iBensbv et al.ll2007h . 

Similarly, Fig. l20l shows the o-Fe vs. [Fe/H] plane for stars 
belonging to the irfm sample. Stars are separated by their ro- 
tation velocities (V) as depicted in the middle panel showing 
V vs. [Fe/H]. Clearly, this is only a rough criterion for the 
division and this selection is not stringent in targeting single 
disc "components", yet a striking difference appears. At each 
metallicity, the stars with high negative V velocities (open cir- 
cles) have higher average aFe; the difference is indeed small in 
terms of aFe, but statistical significant. This can be expected 
because stars with such a large asymmetric drift should be sig- 
nificantly older than the remaining population (because of the 
asymmetric drift-dispersion and the age-dispersion relations), 
which is confirmed by their age distribution in the lower right 
panel, which is indeed far olde r Our ana l ysis t hus clearly con- 
firms a similar result drawn by iHavwoodl (]2008 ) from a smaller 
spectroscopic sample. Comparison with Fig. [18] also shows the 
correspondence between our identification based on rotational 
velocities and the orbital radii of stars. It is interesting to no- 
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Fig. 19. [Fe/H] vs. aFe for stars in the irfin sample within the 
metallicity calibration ranges and with kinematic information to 
assign statistical membership to the thin or thick disc. Only stars 
with a membership probability higher than 90 percent are shown 
(4655 stars in total). A Gaussian noise of 0.005 dex was added 
on both axes for better displaying all stars. 



tice that lEdvardsson et al.l (1 19931) found a hint that stars with 
high/low orbital radii he on the lower/upper envelope of the 
[ff/Fe] vs. [Fe/H] plot, consistent with what we see here. 

Because of the tight age-metallicity relati on in chemical 
evolution models without radial migration (e.g.. lChiappini et aP 
11997) . older and alpha-richer stars are expected to be more 
metal-poor Yet Fig.|20]rather tells the opposite, with the old and 
alpha-rich stars also being on average more metal-rich than the 
population with high rotation velocities (filled diamonds). The 
emergence of a metal-ri ch, old thick disc was alrea dy present in 
spectroscopic sample of iFeltzing & Bensbvl (12008'). This appar- 
ently surprising behaviour is however readily explained if there 
is no strong age-metallicity dependence, as is the case in ra- 
dial migration models, and if the lagging metal-rich population 
comprises -to some extent at least- objects from the inner disc, 
which are more metal-ri ch thanks to the Galactic metallicity gra- 
dient (ISchonrich & Binnev.,2009a) . 

7. Conclusions 

Low mass, long lived stars are crucial witnesses of the chemical 
and dynamical evolution of the Milky Way, but to properly har- 
vest this information, we must ensure that we have determined 
their astrophysical parameters to the highest accuracy possible, 
given the observational limitations. The Geneva-Copenhagen 
Survey provides the ideal database to achieve this goal: it is 
kinematically unbiased, all its stars have highly homogeneous 
Stromgren photometry, from which stellar abundance informa- 
tion can be readily derived and merging this catalogue with 
Tycho2 and 2MASS provides the multi-band optical and infrared 
photometry needed to derive Teff via the infrared flux method. 

We have carried out a revision of the GCS not only bene- 
fiting from the latest developments in setting the zeropoint of 
the effective temperature scale, but also improving upon the ho- 
mogeneity of the stellar parameters for all stars in the sample. 
In comparison to previous GCS calibrations, our effective tem- 
peratures are hotter; at the same time the improved method- 
ology often reduces the intrinsic uncertainty per star to below 
100 K. This leads to a much better agreement between stars and 
isochrones in the HR diagram, which allows us to directly derive 
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Fig. 20. Panel a): aFe for stars in the irfm sample having Vlsr - V +Vq > 20 km/s (filled diamonds) or < -40 km/s (open circles). 
Only [Fe/H] > -0.73 dex were selected to avoid a metal-poor tail in stars with negative velocities (the cut in [Fe/H] was selected 
to be the same as the lower-most value encountered in stars having Vlsr > 20 km/s, but its exact choice is anyway irrelevant 
for the discussion). A few stars with clear halo kinematics were also excluded. Lower/upper continuous line connects the mean 
aFe in different [Fe/H] intervals for filled diamonds/open circles. Error bars (slightly offset in abscissa for clearer comparison) are 
the standard deviation of the mean in each [Fe/H] bin. Panel b): Vlsr as function of metallicity for all stars in sample irfin with 
kinematic information. Filled diamonds and open circles as in the previous panel. Panel c) and d): normalized metallicity and age 
distributions for the two previous group of stars having Vlsr > 20 km/s (continuous line) or < -40 km/s (dashed line). The value 
Vq - 12.24 km/s was adopted (.Schonrich et al...201O) . 



ages via a Bayesian approach. Because we did not make use of 
metallicity-dependent temperature shifts to reconcile isochrones 
with data, the risk of introducing an artificial age-metallicity re- 
lationship is reduced. Since the adopted effective temperature 
scale has immediate consequences on abundances, we recali- 
brated Stromgren indices versus stellar metallicities using a sam- 
ple of nearly 1500 stars with high-resolution spectroscopic abun- 
dances derived adopting Teff consistent with ours. We thus warn 
that when comparing our results with other studies, it should al- 
ways be kept in mind that differences in metallicities could sim- 
ply reflect the different Teff scales adopted. As a consequence, 
the mean metallicity of previous GCS analyses is increased by 
~ 0.1 dex, now peaking at [M/H] ~ -0.01 dex and thus mak- 
ing the Sun a completely average star given its metallicity (see 
also lAsplund et al...2009.) . It is intriguing to note that in the past 
the higher metallicity of the Sun comp ared to local dwarfs was 
used in support of radial migration (e.g. 'W ielen et al.lll99 6^: our 
analysis suggests that the Sun is not atypical, at least in metallic- 
ity. Instead, we derive other atypical properties for disc stars, to 
explain which, radial migration could be a relevant ingredient. 

For the first time we are able to derive [ff/Fe] estimates from 
Stromgren photometry (named ffFe for the sake of clarity). The 
method becomes more unreliable for increasingly hotter objects 
and also for metal-poor stars (roughly below -1 dex), but fo- 
cusing on disc stars gives a reasonable guidance on the relative 
alpha enhancements for the whole sample. The ability to reach 
this tentative distinction enabled us to bring the metallicity cal- 
ibration to significantly better accuracy by reducing the uncer- 
tainty in [a/Fe] enhancement. The new metallicity scale was 
then checked against open clusters and a moving group, showing 
indeed a high degree of internal consistency with a suggested in- 
trinsic scatter below 0.10 dex in [Fe/H]. The recently measured 
uvby solar colours finally corroborate the agreement between the 
temperature and metallicity scales. 



Having this at hand, we revised and complemented the 
largest existing sample of F and G dwarf stars in the solar neigh- 
bourhood that is kinematically unbiased and gives information 
on ages, the abundance plane, and kinematics. A preliminary 
analysis of this dataset supports a scenario with a strong inter- 
play among those three characters: the metallicity distribution 
function shows increasing broadening at older ages, suggest- 
ing that its wings could mostly comprise stars bom at various 
Galactocentric radii and migrated at the current position over dif- 
ferent timescales. This scenario could also account for the radial 
gradient getting flatter for older ages, though this detection is yet 
uncertain partly because of the short distance baseline covered 
by the GCS, and partly because of the difficulties in disentan- 
gling metallicities, ages, and kinematic selection in the sample. 

A more robust and striking feature comes instead from the 
division of stars in the rotation velocity plane, which are shown 
to have different patterns in the abundance plane and in ages, a 
feature which is unexpected in classical chemical evolution mod- 
els, but seen in spectroscopic studies and naturally explained 
if stellar radial migrations is taken into account. Despite that 
our data show clear support for the radial migration scenario, 
many different processes enter the picture of galaxy formation 
and evolution; future larger surveys will thus be invaluable to 
further constrain the interplay of various scenarios. The results 
presented here are thus an example of the importance of having 
at the same time kinematic, metallicity, and age information to 
uncover the past of our Galaxy. 
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Appendix A: Bayesian age determination 

For the effective temperatures we assume a Gaussian error, 
which is derived for each star as described in Section 12.1.21 
or 12.1.31 for the irfm and clbr sample, respectively. Things get 
slightly more complicated for the magnitude errors. The mag- 
nitude is estimated from the photometric measurements and the 
parallax of these stars, while the latter measurement completely 
dominates the error So, assuming a Gaussian distribution in the 
parallax, we can write 

J -(P-Pq)- 

Pp(p\PQ,crp)^^^e '-i , (A.l) 

where po is the best estimate for the parallax, cto the adopted 
parallax error Converting to magnitude space we thus have 

Pv(p(V)\Vo,crp) ^ Pp(V(p)\V^,,ap)^ = 

where AV - V - Vo is the difference between the magnitude 
V and the best parallax-based estimate for the magnitude Vo, 
o'pr - o-pl Po is the relative parallax error and k is some normal- 
ization constant. For small relative parallax errors this treatment 
does not imply any significant changes, because the error distri- 
bution approaches a Gaussian. However, for ctj,,- > 0.1, the V 
magnitude distribution becomes increasingly skewed, because 
the lower parallaxes produce an extended ta il towards brighter 
magnitudes (see also lCasagrande et al.ll2007l for an analytic es- 
timate of the bias in case of low parallax errors). 



It is crucial to correct the estimated metallicities for system- 
atic biases. In the wings of the MDF more stars are scattered out 
from the more densely populated central regions, so that there 
is a net bias to be expected in the metallicities in the wings 
of the distribution, e.g. stars on the high-metallicity wing have 
on average overestimated metallicities, while stars on the left 
wing of the MDF are expected to have an increased fraction of 
metallicity underestimates. Because of these shifts, a naive use 
of the measured metallicities would introduce an age underesti- 
mate on the high-metallicity side and an age overestimate on the 
low-metallicity wing, which would give rise to an artificial age- 
metallicity correlation. This bias can be reduced and in the best 
case removed by an appropriate metallicity prior that reflects the 
underlying "real" metallicity distribution. It might be tempting 
to use the metallicity distribution itself as metallicity prior in an 
iterative process. However, this is intrinsically unstable, because 
stars would assemble in peaks, growing by attraction of more 
objects. So we took an analytical function that approximates the 
sample distribution: 

/([M/H]) = 

( 387.8m([M/H]) for [M/H] > 0.04 

\ 387.0OT([M/H])cor([M/H]) + 0.8 for [M/H] < 0.04 

with 

/ ([M/H] -0.04)2 \ 
^([M/H])^exp(- ^^ 2.0.122 ) 

cormm = 1 + o.3(e-20([M/H]+o.26) _ ^-6.0^^ 

On the right hand side we simply choose a Gaussian term as 
prior. On the left hand side this is considerably flattened by 
adding the "correction" term. This function has to be multiplied 
with the Gaussian error term in metallicity. Mainly its relative 
slope decides about shifts in the adopted probability distribution 
in metallicity. So with the correction term that flattens the distri- 
bution at low metallicities we can hope to reproduce the actual 
data sufficiently well. 

Altogether the probability distributions in each parameter of 
a star are gained by running over the isochrone points. The prob- 
ability distribution in one parameter x, is gained by 

■PviVis\V,crp)G{log(T,n) - log(T,ns.),(rT) 
■/([M/H],.,)G([M/H] - [M/H],.„cr[M/H]) 

■dm,„,isdT.isd[Mlii],is , 

where the sum runs over all isochrone points is. G{{y - yo), cr^) 
is a Gaussian function with y -yoin the counter of the exponent 
and with dispersion cr,., ot,„_,j is the initial mass of the star on 
an isochrone point, ^ is the age prior, d,„, and t/[M/H] denote 
the effective volume covered by an isochrone point, i.e. half the 
distances to its neighbours in initial mass, age and [M/H], and 
subscript is denotes the values of the isochrone point. From the 
probability distribution the maximum likelihood, median, expec- 
tation and 5, 16, 84, 95 percent values of the age are derived. The 
error cr for any given age is defined as half difference of the 84 
minus 16 percent value and the relative uncertainty as the ratio 
of cr over expectation age. Throughout the paper, BASTI expec- 
tation ages and masses are used in the analysis if not otherwise 
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some technicalities used in GCSII, it is difficult to explain all 
trends that arise in the comparison. Some of the breaks (e.g. the 
depletion of stars around 5 Gyr) could arise because isochrones 
in GCSII were shifted to agree with the data. Similarly, we notice 
that stars in GCSII with undetermined lower confidence limits 
are preferentially assigned young ages, which are responsible for 
some of the horizontal stripes seen in the upperpanel of Fig. lA!2] 



Fig.A.l. Comparison between ages derived assuming the 
Salpeter vs. flat IMF. Red dots are stars belonging to the irfin 
sample. 
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Fig.A.2. Comparison between our revised ages on those in 
GCSII for the full catalogue. AAge is in the sense ours minus 
GCSII. Contour levels in the lower panel same as in Fig. |2] 



specified. In all instances a Salpeter IMF was used, and the effect 
of this choice on f (x,) is minimal. Fig. lA.ll shows the difference 
in BASTI expectation ages when a flat IMF is used instead. 

A comparison between our BASTI expectation ages and 
those derived in GCSII is shown in Fig. lA.2l Age determinations 
are subject to many subtleties: because we do not have access to 



